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ABSTRACT 

We compute the expected angular power spectrum of the cosmic Far Infrared Background (FIRB). 
We find that the signal due to source correlations dominates the shot-noise for £ < 1000 and results in 

anisotropies with rms amplitudes (\J 1(1 + l)CV/27r) between 5% and 10% of the mean for I > 150. The 
angular power spectrum depends on several unknown quantities, such as the UV flux density evolution, 
optical properties of the dust, biasing of the sources of the FIRB, and cosmological parameters. However, 
when we require our models to reproduce the observed DC level of the FIRB, we find that the anisotropy 
is at least a few percent in all cases. This anisotropy is detectable with proposed instruments, and its 
measurement will provide strong constraints on models of galaxy evolution and large-scale structure at 
redshifts up to at least z ~ 5. 

Subject headings: cosmology: theory - cosmology: observation - cosmology: far infrared background - 
cosmic microwave background - galaxies: formation - galaxies: evolution 



1. INTRODUCTION 

The recent discovery of the cosmic Far Infrared Back- 
ground (FIRB, |Puget et al. 1996|. JFixsen et al. 1998|,|Dwek| 
et al. 1998| , jSchlegel et al. 1998] , |Lagachc ct al. 1999| ) and 
the determination of its spectrum creates the opportunity 
for a new pro be of structure formation in the high -redshif t 
universe (e.g. Guidcrdoni ct al. 1998| , Blain ct al. 1999a| ). 
One way of improving our understanding of this back- 
ground is to observe the sky with deep exposures at high 
angular resol ution, to determine t he number counts of dis- 



aim of the present Letter is to predict the lowest-order 
statistical property of the FIRB, the two-point angular 
correlation function, using simplified models for the ori- 
gin of the background flux. Our results demonstrate that 
under conservative assumptions, the correlations are mea- 
surable with Planck 1 and other future instruments, such 
as the balloon-borne Far Infrared Background Anisotropy 
Telescope (FIRB AT) proposed by the TopHat 2 group and 
the Bolometric Larg e Aperture Sub-mm Explorer (BLASE, 
Dragovan, M. 1999[). In addition, the signal might be de- 



crete sources ( Hughes ct al. 199^ Barger et al. 1998 , 
et al. 1999J, gmail ct al. 1997], |Holland ct al. 1998T 



Ealc: 



Blair 



et al. 1999b| , pargcr et al. 1999| , [Pugct ct al. 1999| ), 



and t he redshift distribution of the contributin g sources 
(e.g. [Barger ct al. 1999b|, |Blain ct al. 1999cj). Obser- 



vations with the Sub millimetre Common User Bolometer 



Array (SCUBA) camera (Holland et al. 1998) have indeed 
identified point sources at 450/im and 850/im, which ac- 
count for a large fraction ( <; 25%) of the FIRB at these 



tectable in the highest frequency channel of the existing 
BOOMERanG 3 data. 

To calculate the correlation function, we must first un- 
derstand the nature of the source(s) of the background. 
A compelling explanation is thermal emission by the in- 
terstellar dust of high-redshift galaxies, heated by their 
internal optical and ultraviolet (UV) star-light. A hand- 
ful (~ 30) of sources with the expected properties have 
been resolved in several SCUBA images at 850/im, r e- 



rolcngtha. However, a conclusive identification of opti — 
oal oountorparto, or tho dotormination of tho rodohifto of 



ported in several papers (Hughes et al. 1998 



al. 1998, Eales et al. 1999, 3mail et al. 1997, Hollanc 



Barger et 



et al. 1998| , |Barger et al. 1999J , plain et al. 1999b| ) anc 
summarized by Scott & White (1999), and another hand- 
ful at 175/mi with the ISOPHOT i nstrument on boar d the 
Infrared Space Observatory (ISO) ( [Puget et al. 1999[ ). Ac- 
cording to Hughes et al. (1998) and Barger et al. (1999), 
sources brighter than 2m Jy comprise at least 25% of the 
background at 850/im (350 GHz). The observed source 
counts can be parame terized empirically b y a softened 
power-law form (e.g. [Barger ct al. 1999 ). Using this 
form to extrapolate towards the faint end of the source 
count distribution (down to ~ 0.5mJy), one can account 



those sources are still made difficult by tho lack of detailed 



spectral information, and the uncertainty in the SCUBA 
determined positions. 

A complimentary approach is to observe the background 
at relatively low angular resolution in the "confusion limit" 
with the aim of studying its statistical properties. Such 
analyses have proven fruitful in the context of the Cosmic 
Micr owave Backgr ound (CMB), as we ll as for the opti- 



cal (Vogeley 1997) and near-infrared (Jimenez & Kash 
linsky 1997, Kashlinsky et al. 1997) backgrounds. The 



1 Planck: http://astro.estec.esa.nl/SA- general/Projects/Planck/ 
1 TopHat: rittp://topwcb,gsfc,nasa,g< 



3 BOOMERanG: 



http://www.physics.ucsb.edu/boomerang 
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for the entire background. These conclusions are further 
supported by the somewhat less reliable detections of faint 
sources down to ^TmJy. One should bear in mind, how- 
ever, that such extrapolations are highly sensitive to the 
properties of the high-z population. 

The sources identified by SCUBA are possibl y galaxies. 
Two have been conclusively identified as such ( Frayer et 
al. 1S|98|). Hughes et al. (1998), observing in the Hubble 
Deep Field (HDF), identified tentative optical star-burst 
galaxy counterparts for all five of their sources brighter 
than 2mJy; four of them with redshifts between 3 and 
4 and one with z ~ 0.9. These identifications are not yet 
conclusive, due to the coarse angular resolution of SCUBA, 
and the possibility remains that many of the sources are 
dust-enshroud ed quasars, rather than galaxies (see, e.g. 
Sanders 1999] ). 

It is also not yet clear that the FIRB is entirely due to 
the type of sources seen with SCUBA. Other possibly sig- 
nificant contributors include emission from inter-galactic 
dust, ejected from galaxies by radiation pressure or s uper- 
nova winds (Aguirre 1999, Aguirre fc Haiman 1999[ ), nu- 
merous low surface brightness, dusty protogalaxies, and, 
more speculatively, radiative ly-decaying massive particles 
(Bond, Carr & Hogan 1986). These contributions to the 
background will not be seen by experiments that are tar- 
geting relatively bright point sources. They would all, 
however, contribute not only to the mean level of the back- 
ground, but also to its fluctuations on large angular scales, 
as long as their infrared emission traces the spatial mass 
fluctuations to some degree. 

In this paper, we do not model the galaxy evolution 
process in any detail, but instead adopt a toy model that 
allows us to explore the dependence of the FIRB angu- 
lar power spectrum on the nature of the sources, their 
redshift distribution, and cosmological parameters. More 
sophisticated models exist in the literature that general- 
ize semi-analytical galaxy evolution schemes to make de- 
tailed predi c tions in the infrared regime (e.g. roffolattj 



discrete source detections discussed above, as well as de- 
terminations of the redshift- evolution of the global aver- 
age star formation rat e (e.g. |Madau 1999] [M99]), the UV 
background at z = ([Bernstein 1997), and at 2 <, z <y 4 



(Giallongo et al. 1996[) now allow a more focused discus- 



et al. | 1998[ |Guidcrdoni et al. 1998] JBlain et al. 1999a| ) 



However, even these models are very much driven by the 
data and the semi-analytic approach has resulted in a wide 
range of predictions in the past. Here we choose a sim- 
plified approach, which is sufficient to illustrate the de- 
tectability of the clustering signals, their dependence on 
the large-scale distribution of matter and aspects of galaxy 
formation, and to demonstrate the need for further work. 

An approach even simpler than our own is that of Scott 
and White (1996) who estimated the angular correlation 
function of the FIRB at 850/ito by assuming that the FIRB 
sources have an angular correlation function lik e that of 
Lyman-break galaxies (Giavalisco et al. 199S). They 
viewed the FIRB as a contaminant of the CMB data, and 
were interested in a rough estimate of what these corre- 
lations might possibly be. As a result, they could ignore 
the fact that the projection effects would be different for 
the two different classes of objects. In the present paper, 
we investigate the FIRB as an interesting signal in its own 
right, and model it in sufficient detail to understand the 
dependence on large-scale structure and galaxy formation. 

The spectra and correlation functions of both the near 
and far infrared backgrounds have been considered in gen- 
eral terms in the pioneering works of Bond, Carr & Hogan 
(1986, 1991). The recent measurements of the FIRB, the 



sion. When we calibrate our models using the recent in- 
frared data, we find that the clustering results in contrasts 
of about 10% in the FIRB, which is sufficiently strong 
to dominate the shot-noise (estimated from the SCUBA 
detect ions) and to be detectable b y proposed future mis- 
sions, limenez fc Kashlinsky 1997| have calculated the an- 
gular power spectra for the near infrared background, and 
also find contrasts of ro ughly 10%. Indeed, very recently, 
Kashlinsky et al. 1999 have claimed a tentative detection 
of fluctuations of cosmic origin in the DIRBE (Diffuse In- 
frared Background Explorer) maps from 1 to 5 microns, 
with the expected amplitude, or possibly slightly higher. 
The angular power spectrum in Scott and White (1996) is 
bracketed by those in our range of models, although most 
of our models tend to have somewhat higher amplitudes. 

Our work has been strongly motivated by future mis- 
sions including Planck and the FIRB AT. Both are capable 
of detecting these fluctuations, the FIRB AT by observing 
in eight channels with central wavelengths ranging from 
230/im to 940/im with per-channel sensitivities between 
25 and 130 /j,K s 1 / 2 , and an angular resolution of 6'. The 
broad range of frequencies is important for separating the 
FIRB fluctuations from those of dust in our own galaxy. 
It is possible that this signal is discernible in the existing 
BOOMERanG data, although the limited spectral cover- 
age will make discrimination from dust in our own galaxy 
difficult. 

The rest of this Letter is organized as follows. In § 2, 
we describe our toy model for the mean FIRB. In § 3, we 
extend this model to include the fluctuations, by assum- 
ing cold dark matter (CDM) power spectra, and that the 
FIRB light is a (biased) tracer of mass. In § 4, we show 
results in several models, including those with a uniform 
bias, and a redshift-dependent bi as calculated accor ding 
to a prescription for galactic halos (Mo & White 1996). In 
§ 5, we discuss the results, and what can be learned from 
planned observations of the FIRB on large (> 5') angular 
scales, and finally in § 6, we summarize our conclusions. 

2. MODELING THE MEAN FIRB 

In our model, the FIRB arises from thermal dust emis- 
sion. To compute the mean level of the flux, the main 
ingredients are the mass density f2d, and temperature T^ 
of dust, and the evolution of these quantities with redshift. 
The angular fluctuations depend further on the spatial dis- 
tribution of the dust, as we will discuss in § 3 below. In 
order to compute the evolution of Cla(z) and T<i(z), we 
rely on the recent determinations of the evol ution of the 
globa l average star-formation rate (SFR, see [Lilly et al 



1996, M99 and references therein). In particular, we as- 
sume that both the UV emissivity that determines the dust 
temperature, and the rate of dust production, are propor- 
tional to the global SFR (corrected for dust absorption), 
p*(z), as given by M99. Our motivation for relying on the 
star-formation history is that the stellar UV flux is known 
to dominate that of quasars at all redshifts z <, 5, with 
quasar contributions less than ~ 20% (M99). We assume 
further that the dust has a composition of 50% graphites 
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and 50% silicates by mass, and compo site cross-sections as 
in the Milky Way (Draine & Lee 1984). Finally, we assume 
that the UV emissivity has the spectral shape obtained by 
summing individual zero-age main sequence s tellar spec- 
tra w eighted by a Scalo (1986) mass function (Bruzual & 
Chari ot 1996J ). 
W 

Table 1 

The parameter values in our standard model and its variants. The 

last two columns are the values of the dust density and 

temperatures at z = (in Kelvin), determined from fitting the 

FIRB at 850/tm. See text for discussion. 



Under these assumptions, the dust temperature is de- 
termined by the UV radiation, i.e. by requiring that the 
amount of energy absorbed in the UV is equal to that ther- 
mally radiated in the infrared. The proportionality con- 
stants for the dust production rate and the UV emissivity 
are then determined by fitting the observed spectrum of 
the mean FIRB. The relevant equations ar e summarized in 
§ 3 of Aguirre & Haiman (1999) (see also Loeb fc Haiman| 
1997|] Note that we have not assumed that the dust is per- 
fectly uniformly distributed, but only that wherever there 
is dust, it is exposed to the same radiation field. Also note 
that in this scenario, as long as we are interested only in 
the mean level of the FIRB, we do not need to specify 
the spatial distribution of the dust, or the type of source 
from which the infrared emissivity arises. To get the an- 
gular power spectrum, we will assume below that dust is 
a (biased) tracer of mass on large scales (see § 3). 

The parameters of our models are summarized in Ta- 
ble ||. Our standard model is a ACDM cosmology, with 
the SFR taken from M99. In this model, we find that in 
order to fit the FIRB, the graphite/silicate dust tempera- 
tures at z = need to reach 14. 4K and 10. IK, while the 
total mass density of dust has to reach £1h = 2.7 x 10 -5 . 
The redshift evolution of these quantities are shown in the 
upper panel of Figure [jj In particular, Qd rises contin- 
uously as dust is accumulated from the increasing num- 
ber of stars. The dust temperatures deviate from the 
CMB at z « 7, and stay roughly constant (silicates at 
~15K, graphites at ~25K) until z ~ 1, at which point 
the SFR drops sharply, the dust heating becomes less ef- 
ficient, and the dust temperatures start dropping again. 
The dust temperatures we obtain in our fiducial model 
at z — are somewhat colder than, e.g. the tempera- 



ture of dust in the Milky Way (~ 16 - 20K, [Reach et 
al. l |995| ) . The temperatures are higher in our "flat" 
SFR model (~ 17/13K, see Table p] and discussion be- 
low), and are in better agreement with the Milky Way 
value. We therefore note that our conclusions are not 
dependent on an accurate reproduction of dust temper- 
atures in local galaxies. In the lower panel of Figure [j], we 
show the redshift evolution of the comoving dust emissiv- 



ity jd(z) = Pd x «d,j/(i+2) x B v(i+z) PV(1 + z)] in units 
of 10 -21 erg cm -3 s _1 sr _1 Hz -1 at the (observed) wave- 
lengths of 850/xm (lower solid curve) and 450/zm (lower 
dashed curve). Here k is the dust opacity coefficient and 
B v is the Planck function. The curves in the figure il- 
lustrate the faster fall-off towards higher redshifts of the 
emissivity at 450/im compared to 850/im, due to the fact 
that the peak of the grey body dust observed at 450/im 
occurs at lower z than for the longer wavelength 850/im. 
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Fig. 1. — Upper panel shows the evolution of dust temperatures 
in our standard model, compared with the CMB temperature, as 
well as the mass density of dust. The labels on the right axis refer 
to logf2 d . Lower panel shows the comoving dust emissivity j' d at 
450/im and 850/im (bottom curves). The top curves show the quan- 
tity 6j' d , where b is the bias parameter of galaxy sized (10 12 Mq) 
halos. 

We considered several variations on our standard model, 
summarized in Table [lj. To explore the dependence of our 
results on the underlying cosmology, we have evaluated an 
open CDM model with Sl m = 0.3, and a standard CDM 
model. We then studied a "flat" SFR model, in which the 
SFR is constant in redshift. The motivation for consid- 
ering this model is that the UV flux that heats the dust 
inside an individual galaxy may not evolve in the same 
way as the global average SFR does. Indeed, the evolu- 
tion of the global SFR is driven mainly by the change in 
the number of galaxies, r ather than th e evolution of the 
characteristic brightness ( Madau 1997 ). In this case, the 
dust temperature within each galaxy might stay roughly 
constant, provided it is heated by long-lived stars. To 
mimic this scenario, we have included a decreasing overall 
dust abundance according to the original M99 SFR, but 
set the UV emissivity to be a constant at all redshifts. 

As an alternative to this scenar io, we also c onsidered a 
"hot" dust model, adopted from Blain et al. 



199jkj . In 

this model, the dust in each galaxy is kept hot by the 
more intense radiation from short-lived, intense bursts of 
star-formation, after which the dust cools down, and its 
emission is subsequently ignored. The total dust emis- 
sion at a given redshift in this case is proportional to the 
UV emissivity, which, in turn, is proportional to the SFR. 
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Following Blain et al. 1999c , we assume a constant dust 
temperature of 40K, and a star-formation rate that is, in 
effect, a sum of the M99 SFR, and an additional Gaus- 
sian peak centered on redshift z = 2.1. Note that in this 
model, only a fraction of the total dust abundance is "lu- 
minous" at any given time - the entry for fi^ in Table [0 
refers to the maximum luminous dust density, reached at 
z = 2.1. Finally, we consider a "high- z" variation on our 
standard model, in which we postulate an additional peak 
in the SFR at redshift z = 7, with SFR oc exp[-(z-7) 2 /2], 
and an amplitude such that this hypothetical high-redshift 
population accounts for 50% of the FIRB at 850/zm. The 
choice of this large burst of star formation is somewhat 
ad-hoc. Our motivation for considering this scenario is to 
characterize the possibility that the FIRB has a significant 
contribution from very high redshifts - an option that can 
not be ruled out by present observations. Note that sig- 
nificant star-formation at redshifts z <; 5 can still have 
escaped detection by, e.g., HST. 
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Fig. 2. — The spectrum of the FIRB predicted in our standard 
model (solid curve). The dot— dashed curves show the measurements 
with ±1(T uncertainties from Fixsen et al. (1998), and the long 
dashed curves show the CMB for reference. The data— points are 
taken from Schlegel et al. 1998. 

In Figure @ we show the spectrum of the FIRB in our 
standard model. The dot-dashed curves show the fit to 
the measurements by Fixsen et al. (1998), while the solid 
line shows our model prediction. The long-dashed lines 
show the CMB, and the vertical solid lines mark the ob- 
servational wavelengths of 450^m and 850/im. As the fig- 
ure reveals, the prediction from our model is in reason- 
ably good agreement with the FIRB; in particular, it is 
within the quoted ±ltr error bars. The quality of fits in 
the other models from Table |l| is similar to that of the 
standard model shown in Figure H. The real sources of the 
FIRB are much more complex than our simplified model 
assumes. The UV flux density and the dust density are 
not expected to be spatially uniform, and dust in differ- 
ent galaxies could have different temperatures, broadening 
our predicted dust-emission peak. However, our model is 
sufficient to allow for a rough prediction for the angular 



power spectra at various wavelengths, and to provide the 
framework for a discussion of what might happen with 
more realistic models. 

3. FIRB FLUCTUATIONS 

The model described so far accounts for the amplitude 
of the FIRB, but makes no specific reference to the spatial 
distribution of the infrared-emitting dust. In this section, 
we characterize the spatial distribution in order to derive 
the angular power spectrum of the FIRB. One simple as- 
sumption is that "dust traces mass", i.e. that the local 
density of dust is proportional to the total mass density 
at every point in the universe. In reality, most of the 
dust could be confined in galactic halos; in which case the 
spatial distribution of dust emission would follow that of 
galactic halos. The correla tion function of dark halos is 
related ( Mo fc White 19"96| ) to that of mass by the bias 
parameter b = b{M\ la _\ , z). Here we adopt the modified 
formula of Jing (1999), which has been shown to repro- 
duce the results of numerical simulations on a wider range 
of scales (including -Mhaio < M*). Under these assump- 
tions, the bias b is scale-independent for a fixed halo size, 
but evolves in redshift. The correlation function between 
pixels separated by angle 9 and at frequencies v and v' is 
given by: 



0/ i 1 

C vu (6) = J2 =^jz-Cr Pi (cos 6>) 



i 



4tt 



Of = 1 f k 3 P(k)f ul (k)f ull (k)f 



(1) 
(2) 



Mk)= ji(kr)b(r)D{r)j d (v,r)a(r)dr (3) 



where r = c(t]q — rj) is the coordinate distance of an event 
at conformal time r\ on our past light cone, 770 is the con- 
formal time today, P(k) is the power spectrum of the 
matter today, a(r) is the scale factor normalized so that 
a(0) = 1, D(r) is the linear theory growth factor, and ji is 
the spherical Bcsscl function. In a matter-dominated uni- 
verse, D(r) — l/a. For the more general case of non-zero 
curvature and/or non-zero cosmological constant we use 
the fitting formula of Peebles (1980) and Caroll, Press & 
Turner (1992), respectively. 

Equations 2 and 3 are a version of Limber's equation 
( Limber 1953| , Peebles 1980), although with power spectra 
instead of correlation functions and generalized to describe 
correlations between different components. Note that pix- 
els at unequal frequencies will not be perfectly correlated 
because jd(v, z) and jdiy' , z) are not proportional to each 
other for v ^ v ' . The unequal-frequency correlation func- 
tion (at 8 = 0°, with 5' smoothing) has been calculated by 
Bouchet & Gispert (1999) by using simulated Planck maps 
generated using the semi-analytic model of Guiderdoni et 
al. (1998). 

As equations (1-3) show, Ci scales roughly as b 2 . To 
illustrate the effect of a non-negligible bias, in the lower 
panel of Figure fil we show the evolution of the comoving 
emissivity jd(v,z) at 450^m and 850/im (bottom curves), 
together with the quantity jdiy, z)6(Mh a io, z) (top curves). 
We have assumed Mh a i = 10 12 M Q , which is valid if dust 
emission arises from dark halos of this size (roughly the 
size of galactic halos, as well as the host halos of typical 
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quasars, see Haiman fc Loeb 1998 ). The figure shows that 



biasing significantly boosts the contribution to the signal 
from z <; 1, but has negligible effect at Z ^ 1. In what 
follows, we alternately use this prescription for biasing and 
a time-independent bias of b — 3 , which is roughly that of 
Lyman-break galaxies at z ~ 3 (Giavalisco et al. 1998). 
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Fig. 3. — Predicted angular power spectra in units of (/iK) 2 (an- 
tenna temperature) for different cosmologies: ACDM (solid), SCDM 
(dashed) and OCDM (dotted), all assuming M99 SFR. The upper 
panels assume linear bias with a constant 6 = 3, and the lower pan- 
els assume a z— dependent bias. The left panels show the results 
at 450/im and the right panels at 850/im. The long-dashed line 
on the upper left panel shows the sensitivity of FIRBAT, assuming 
bins of 51 = 50 and a one week flight. For reference, the CMB power 
spectrum is also shown (dot-dashed lines in the 450/im panels). The 
light solid lines (with bounding dashes for the 450/im panels) are es- 
timates of the shot-noise power with indications of the uncertainty. 
We do not attempt to characterize the uncertainty at A = 450/tm. 

For our standard (ACDM) model, we have adopted the 
fitting formulae for the power spectrum given by Eisenstein 
& Hu (1999) with power-law index n = 1 and normaliza- 
tion CT8 = 1, while in the Open and SCDM models, we used 
(n = 1.3, as = 0.85) and (n = 0.7, a 8 = 0.6), respectively. 
These values for n and as were c hosen to roughly agree 
with both the clus ter abundance (Viana fc Liddlc 1999 ) 
and COBE-DMR ( |Bunn fc White 1996Q constraints. For 
the conversion between coordinate distance and scale fac- 
tor for the ACDM model, we used the fitting formula of 
Pen (1999). 

The total C/ is the sum of that due to correlations, which 
we have calculated, and that due to the discrete nature 
of the sources. This "shot-noise" Ci was calculated by 
Scott & White (1999) at A = 850/um who used the double- 
power law LF constrained by SCUBA data mentioned ear- 
lier and is shown in Fig. 3 rising like l 2 (as all white noise 
power spectra do). The LF at A = 450^m is much more 
uncertain; we estimate the shot-noise at this wavelength 
simply by assuming that the shape of the LF is indepen- 
dent of frequency. Since the antenna temperature of the 
mean FIRB is roughly the same at these two wavelengths 
(T mean ~ 60/iK) , this means that the shot-noise at 450/im 
is the same as at 850/mr. 



4. RESULTS 

Our results are plotted in Figures [| and [|. Figure || 
shows the predicted angular power spectra in units of 
(/iK) 2 (antenna temperature) for different cosmologies: 
ACDM (solid), SCDM (dashed) and OCDM (dotted), all 
assuming M99 SFR. The upper show our models with a 
constant linear bias 6 = 3, and the lower panels show mod- 
els with a z-dependent bias. The left panels show the re- 
sults at 450/im and the right panels at 850/im. The long- 
dashed line on the upper left panel shows the sensitivity of 
FIRB AT 1 assuming bins of 81 = 50 and a flight duration of 
6 x 10 5 s. For reference, the CMB power spectrum is also 
shown (dot-dashed lines) in the 850/im panels. The light 
solid lines with bounding dashes in the same panels are 
estimates of the shot-noise power with indications of the 
uncertainty. We do not attempt to characterize the un- 
certainty at A = 450/im. Figure |] demonstrates how our 
results change due to variations away from our standard 
model (solid line): we show the flat SFR model (dotted), 
hot dust model (dashed), and high- 2: peak model (dot- 
dashed). 



850/im 



450/im 




10 100 1000 



100 1000 



Fig. 4. — Angular power spectra aslin Figure H, but showing vari- 
ations away from our standard model (solid line): flat SFR (dotted), 
hot dust (dashed), and high-z peak in SFR (dot-dashed). 



The 



first thing to note 



is that 
150 and : 



the contrast 
1000 is about 



l)Cj/(27r)) between £ 
10% of the 6O/1K mean and is mostly due to the intrinsic 
clustering of the FIRB, rather than shot-noise. We can 
roughly understand this basic result analytically with the 
following expression: 



i(i + m _ ^ 



2tt 



rpl 



(1 + Zpcak) 



A 2 (fc) 



ir/k 



(4) 



where A 2 (fc) = k 3 P(k)/(2ir 2 ) is the contribution to the 
variance of the matter density from each logarithmic in- 
terval in fc, a is the coordinate distance over which there 
is substantial emission, and fc ~ V r peak where r pea k and 
Zpoak are the coordinate distance and redshift of the peak 
emission. The right-most term is the inverse square of the 
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"wash-out factor" arising from the incoherent summing of 
structure at different redshifts, which reduces the contrast. 
Taking z pca k — 1 and I = 1100 we find k ~ O.UhMpc^ 1 and 
A 2 (k) = 3 for our standard model. Further taking 6 = 3 

and a ~ 3500 h^Mpc we find [ ' ( ^ C ' ]i = noo ^ 40//K 2 . 
The shape of 1(1 + l)Ci/(2ir) (rising from low I, plateau- 
ing, and then slowly dropping) is due to a) the fact that 
A 2 (k) oc k 3+n at low k, flattening out to k n ~ x at high k 
and b) the wash-out factor resulting in another factor of 
fc- 1 . 

Note that the only quantities that matter for Ci 
are the matter power spectrum, a(r), and the product 
b(r)D(r)jd(r). With the matter power spectrum and a(r) 
fixed, any effect that pushes the bDj product to peak at 
higher redshifts, shifts the spectrum to higher I. Thus, for 
the constant bias cases, the flat SFR model has slightly 
more power at low I and the high-z SFR model is shifted 
to higher I. Likewise, the hot dust model has the steep- 
est rise from low I to high I because its emissivity at both 
450/im and 850//m rises from low z more quickly than any 
of the other models. 

The redshift-dependent bias, b(z), increases monoton- 
ically and at z —0, 2 and 7 has values of 1, 3 and 18 
respectively. Because b(z < 2) < 3 the redshift-dependent 
bias generally decreases the fluctuation power compared 
to the b = 3, z=independent ^sc This is pspppially txue 



high resolution and masking out pixels with fluxes above 
some threshold. Fluctuation analyses of point-source 
cleaned SCUBA maps have been done by two groups 
dHughes et al. 1998| , |Bo"rys, Chapman fc Scott 19981 ). How- 
ever, at these sub-arcminute angular scales, the wash-out 
factor is very large, and the correlated component is ex- 
pected to be very small. These analyses were motivated by 
the desire to see the shot-noise due to unresolved sources 
below the flux cut, and resulted in constraints on the faint- 
end slope of the number counts. 

An actual detection of the FIRB fluctuations will be 
highly valuable. The spectrum of the FIRB anisotropy 
may eventually be better known than the spectrum of the 
FIRB mean due to the experimental advantages of differ- 
ential measurements. Such an improved determination of 
the FIRB spectrum would provide detailed constraints on 
galaxy formation models, in addition to those from the 
amplitude and scale-dependence of the anisotropy signal. 
It will also be interesting to examine the cross-correlations 
of the FIRB fluctuations with other data such as the 
CMB (FIRB sources may lens the CMB, or both back- 
grounds may be lensed by the same mass distributions), 
radio sources (from, e.g., the FIRST VLA survey or the 
MAP 4 22 GHz channel), galaxies , or quasars at diff erent 
rcdshift slices from the 2dF (e.g., Folkes et al. 1993) and 



Sloan Dig ital Sky Survey (SDSS, see e.g., Gunn fc Wein 



at lo w / and also for A = 450/;,m. The exceptions are the berg 1995[ ), and the X-ray, and near infrared backgrounds 



high-z SFR case and the hot dust case at A 
to the significant emission from z > 2. 

5. DISCUSSION 



850/im, due 



Our simplified model of the sources is insufficient for 
many purposes, e.g., predicting the luminosity function, 
but is adequate for predicting the large-angular scale 
power to within factors of ~ 2. The key assumption is 
that the dust-light is a biased tracer of the mass; i.e., 
where there is more mass there is proportionately more 
light. Note that in principle, a non-negligible fraction of 
the unresolved FIRB at 450/im, and especially at 850/im, 
could ari se from the direct emi ssion from optically faint 
quasars ( Haiman fc Locb 1998] ); our approach would be 
equally applicable in this case. The fluctuation level of 
10% is fairly model- independent, although linearly depen- 
dent on the bias of the sources. If th e bias is unity, as i s 
the case for galactic halos near z « ( Mo fc White 199G| ) , 
then the contrast may be as small as a few per cent. 

Contrasts observed to be smaller than a few percent 
would be surprising. This would indicate that most of the 
FIRB comes from sources with b < 1, so that b(r)D(r) is 
small where j(r) is large. Examples of such sources would 
be galaxies near z ~ 0, or intergalactic dust at high red- 
shifts. The former case might be difficult to reconcile with 
the relatively low number of nearby IRAS sources. In the 
latter case, one would need to invoke dust temperatures 
higher than expected based on the estimates of the UV 
background ( Bernstein 1997| ), in order to get the correct 
spectr um for the mean FIRB in the full 150 — 1000/im 
range (Aguirre & Haiman 199E). 

The correlated component can even be measured on an- 
gular scales where it is smaller than the shot-noise. The 
shot-noise contribution can be reduced by observing at 

4 Mj4P: pttp://map,gsfc, nasa.gov 



(the same sources could contribute to the backgrounds at 
both of these wavelengths). 

Further development of predictions for the statistical 
properties of the FIRB is clearly warranted. The correla- 
tion function, the statistical property we have focused on 
here, may prove easier to understand than the luminosity 
function — the bright end of which is probably dominated 
by very rare events (Bond 1999). Correlation function 
predictions will be refined by detailed modeling that is in- 
formed by additional observations. One useful input will 
be the clustering properties of galaxies and quasars, mea- 
sured accurately in forthcoming redshift surveys, such as 
2dF and SDSS. It may then be possible to infer the na- 
ture of the sources of the FIRB by comparing its statistical 
properties to those of galaxies and quasars. 

In the preceding, we have implicitly assumed that the 
large-scale distribution of matter (P(k) and its evolution 
with redshift) will have been determined by high-precision 
CMB anisotropy measurements and redshift surveys, and 
have therefore been focusing on the dependence of Ci on 
the nature of the sources. However, we emphasize that the 
FIRB Ci is sensitive to large-scale structure at redshifts in- 
termediate to those that will be directly probed by these 
redshift surveys (z < 1) and CMB missions (z ~ 1100). It 
will also be sensitive to wavelengths too small to be con- 
strained by the CMB measurements, and at any given co- 
moving wavelength the matter fluctuations will be better 
approximated by linear theory than they are at lower red- 
shifts. The relation of Ci to large-scale structure is com- 
plicated by its simultaneous dependence on jd( v i z) and 
the bias properties of the sources, but these complications 
will be reduced both by improved theoretical modeling of 
the sources and also by high-resolution, deep observations, 
together with identification of counterparts at other wave- 
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lengths. Thus we view point source observations and mea- 
surements of fluctuations on large angular scales as com- 
plementary: the recovery of the power spectrum on large 
scales from Ci observations is aided by measurement of the 
point sources, which is incapable of determining the power 
spectrum alone. 

6. CONCLUSIONS 

The recently discovered cosmic Far Infrared Background 
and the galaxy counts at 850/im have opened a new 
wavelength at which galaxy formation and evolution and 
large-scale structure can be studied empirically. A signif- 
icant fraction of the FIRB has already been resolved by 
SCUBA into discrete sources; however, due to the lack 
of detailed spectral information, identified counterparts in 
other bands, and secure redshift determinations (except 
for the two sources in Frayer et al. 1998), fundamen- 
tal questions still remain unanswered. Does the popu- 
lation of the observed sources indeed account for the full 
background? What is the nature of the observed sources: 
are they galaxies, dust-enshrouded AGN, or a mixture of 
both? What are their redshift distributions? We have 



argued that measurements of the FIRB correlation func- 
tion can help answer these questions and also that the 
FIRB provides a unique probe of the large-scale distribu- 
tion of matter at intermediate redshifts. We have found 
that under simple, but broad, models for the mean FIRB, 
the anisotropy of the unresolved FIRB intensity is at the 
3 to 10 percent level. These fluctuations are measurable 
with proposed balloon-borne instruments and future space 
missions whose datasets will add new, useful constraints 
on large-scale structure and models for the formation and 
evolution of galaxies and quasars. 
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